Pumping charge with ac magnetic fluxes and the dynamical breakdown of Onsager 

symmetry 
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We study the transport properties of setups with one and two mesoscopic rings threaded by ac 
magnetic fluxes of the form &(t) — & dc + $ ac cos(J7oi + 8) and connected to two different particle 
reservoirs. We analyze the conditions to generate a pumped dc current in the adiabatic regime. We 
also study the symmetry properties of the induced dc current as a function of the static component 
of the flux, $d c , with and without a dc bias voltage applied at the reservoirs. We analyze, in 
particular, the validity of the Onsager-Casimir relations for different configurations of the setups. 
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I. INTRODUCTION 

A mesoscopic ring threaded by a magnetic flux is one 
of the paradigmatic systems to discuss the fundamentals 
of quantum electronic transport. Several experiments^— 
and theoretical works^ devoted to the investigation of 
the Aharonov-Bohm effect and related phenomena origi- 
nated by static magnetic fluxes, constitute milestones in 
this area of Physics. 

In the case of a time-dependent magnetic flux, a very 
interesting example corresponds to that where the flux 
increases linearly in time. In this case a constant elec- 
tric field is induced along the circumference of the ring, 
generating a time-dependent current. This problem was 
introduced by Buttikcr, Imry and Landauer in the early 
times of the theory of quantum transport^ and triggered 
very insightful discussions on the role of inelastic scatter- 
ing as a necessary ingredient to generate a dc component 
in the induced current along the ringer— The discussion 
of these effects, along with the effect of disorder, and 
the comparison of the current generated in this setup 
with the current generated by a dc bias voltage were the 
subject of further investigations *£r— While these studies 
have been illuminating from the conceptual point of view, 
the experimental implementation of a magnetic flux with 
a linear dependence in time is not very practical. Instead, 
fluctuating magnetic fluxes are much more usual in the 
laboratories. 

A magnetic flux with an harmonic dependence on time, 
characterized by a frequency f^o, threading a ring con- 
nected to a single wire was considered by Biittiker in Ref. 
9. That work was developed in the framework of the dis- 
cussion of the combined effect of the electromotive force 
(emf), induced by the magnetic flux, and the inelastic 
scattering in inducing a dc current along the circumfer- 
ence of the ring. The coupling to the wire was introduced 
to provide a concrete mechanism for inelastic scattering 
and decoherencc which leads to an induced current along 
the ring with a finite dc-component. More recently, fur- 
ther details of the induced current along the ring were 
analyzed in the same setup^ Other mechanisms to in- 



troduce decoherence in the ring without coupling to elec- 
tronic reservoirs were also considered^ as well as the 
effect of spin-orbit coupling^ In the case of the ring 
connected to a single reservoir, a time-current is induced 
through the contact between the ring and the reservoir 
as a consequence of the driving. Such current has a zero 
dc-component. However, when additional reservoirs are 
coupled to the ring as in the sketch of Fig. la, the driven 
ring may behave as a quantum pump that induces a cur- 
rent with a net dc component between the reservoirs. 
The aim of the present work is, precisely, to analyze the 
behavior of such a current. 

In the last years, quantum pumps have received a 
lot of attention from both experimental and theoret- 
ical communities i2ir— The basic idea of these setups 
is the generation of a dc current in the absence of 
an explicit dc bias. Most of the studied devices are 
mesoscopic structures locally driven by ac gate voltages. 
The key to induce a dc current under these conditions 
is the breaking of time-inversion and spacial-inversion 
symmetries! 25 ' 28 " — However, depending on the mecha- 
nism employed, different behaviors are expected for the 
induced dc current, as a function of Q,q. In particular, 
the so called adiabatic regime, characterized by a dc cur- 
rent with a linear dependence in Oo, is achieved when 
pumping is induced by a setup with two time-dependent 
parameters* 2 ^ which are usually two ac potentials applied 
at different places of the structure and oscillating with a 
phase-lag. Pumping mechanisms have been studied in 
rings threaded by a static magnetic flux and driven by 
applying a local ac gate voltage at some point of the 
circumference i2£r— Adiabatic pumping was proposed to 
be generated by two rings threaded by harmonic mag- 
netic fluxes oscillating with a phase-lag and connected 
by a tunneling contact 41 In the present work we study 
configurations containing one and two ring connected to 
two different particle reservoirs and threaded by a mag- 
netic fluxes with dc and ac components oscillating with a 
frequency Qq. As mentioned before, this system behaves 
as a quantum pump, which induces a dc current between 
the two reservoirs. Our first goal is to identify under 
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which conditions an adiabatic regime is expected in this 
setup. We also extend our analysis to study the trans- 
port properties when an additional small bias voltage V 
is applied at the reservoirs. 

An important related issue is the behavior of the in- 
duced dc current as a function of the static component of 
the flux, $ dc and the corresponding analysis of the valid- 
ity of Onsager-Casimir relations in this setup. The latter 
are consequences of microreversibility and enforce the lin- 
ear stationary conductance of a two-terminal system to 
be an even function of $ rfc . For large voltages V, beyond 
the linear response regime, there is no reason to expect 
that symmetry in the induced current and its breakdown 
has been suggested to have interesting consequencies on 
thermoelectric effects^ Several recent works have been 
devoted to the study of mechanisms for breaking Onsager 
symmetry in the non-linear conductance theoretically 36 
as well as in several experimental settings . 37 ' 38 In most 
of these cases, the source for the asymmetric behavior of 
the current as a function of $ dc was identified to be the 
effective voltage profile induced along the biased struc- 
ture as a consequence of the Coulomb interaction. In 
the case of rings threaded by dc fluxes while biased by 
ac voltages there are also experimental results, which are 
supported by semiclassical theoretical arguments, indi- 
cating that the Onsager-Casimir relations are in general 
not valid for the conductance associated to the rectified 
current. 38 In this context, the second goal of the present 
work is to check the validity of Onsager symmetry in the 
context of setups containing rings threaded by magnetic 
fluxes with dc as well as ac components. To this end, we 
define appropriate conductance coefficients to character- 
ize the dc-current and analyze the symmetry properties 
of these coefficients as functions of $ dc . 

The paper is organized as follows. In section II we 
present the model, the theoretical treatment to evaluate 
the dc currents, as well as the definition of the different 
transport coefficients. In section III we discuss the condi- 
tions to have adiabatic pumping. Section IV is devoted 
to analyze the symmetry properties of the pumped dc 
current as a function of the dc magnetic flux. This anal- 
ysis is extended in Section V, where we also consider the 
effect of a dc bias voltage. In section VI we close with a 
summary and the conclusions. 



II. THEORETICAL APPROACH 
A. Model 

The system we consider is sketched in Fig. [TJ It 
consists in one or two single-channel rings of length 
L threaded by harmonically time-dependent magnetic 



a) 



fluxes of the form 3>j(i) 



3 



■ cos(f2 i + Sj), where 



j = 1,2 labels the different rings, which are connected 
to one-dimensional left (L) and right (R) wires that play 




FIG. 1. Sketches of the systems considered in this work. Panel 
a): a single metallic ring threaded by a magnetic flux $(t) 
and connected to two particle reservoirs with chemical poten- 
tials fiL and \xr respectively. Panel b): two connected rings 
threaded by magnetic fluxes <F/(t) = + <1>° C cos(fioi + 
Sj), j — 1, 2. The rings are described by tight-binding chains 
with N sites. The reservoirs are attached at the sites II and 
Ir are those sites at which the leads are attached through 
tunneling contacts Wk L and Wk R - 



the role of reservoirs. The ensuing Hamiltonian is 



H = H r (t) 



E 

a=L,R 



(1) 



where the Hamiltonians for the rings correspond to tight- 
binding models with lattice constant a and N sites each, 
thus L = Na. It reads 



JV 



H r (t) 



H.c 



3=1 

-w, 



1=1 

[4/12,2 + H.c 



(2) 



where the phases <f>j(t) = $>j(t)/L, j = 1,2 are the Pierls 
factors, which account for the magnetic fluxes threading 
the rings in units of the flux quantum he/ e. For simplic- 
ity, we consider spinless electrons and we impose peri- 
odic boundary conditions N + I = I in each ring. The 
last term represents the coupling between the two rings. 
The leads or reservoirs are represented by non-interacting 



Hamiltonians H n 



„t 



Ck 



while the Hamilto- 



nian representing the contact between the leads and the 
rings reads 



Hr, 



a=L,R k a 



c k Cl c + H.c. 



(3) 



where l a are the sites of the rings at which the leads 
are attached. The case of a single ring corresponds to 
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u>2 — w c = in H r (t) and both l a lying on the first ring. 
In what follows, we adopt units where e = c = K = 1. 
We also set to unit the lattice parameter a. 



B. dc current and transport coefficients 

In the most general case, we assume a small voltage V 
applied in the setup, which is represented as difference 
between the chemical potentials of the L and R reser- 
voirs, respectively, /i£ = /i + V and (Ir = (j,. In order 
to analyze the transport properties, we follow the proce- 
dure of Refs. [U, to which we refer the reader for further 
details. The evaluation of the dc current flowing through 
the contact between the reservoir a and the ring to which 
it is attached cast 



E E 

j3=L,R n 



Z7T 



xT a (uj + nft )\Qi a i g(n,u))\ 2 Tp(u)), 



(4) 



where T a (uj) = 2nJ2k a \ w k a | 2( K W _ £ k a ) and f a = 
©(Ma — w) is the Fermi function (we focus on zero tem- 
perature) corresponding to the reservoir a in equilibrium. 
The retarded Green's function for sites 1,1' of the rings 
(for simplicity we use a single label I to identify the site 
and the ring) is expressed in terms of the Floquet-Fourier 
representation 



Gf M (t,t') 



du) 
2^ 



Gft / (i, W )e-*"<*-*'> ) 



^2e- innot g lM (n,u;). (5) 



To evaluate the latter Green function, it is convenient to 
express the Hamiltonian for the rings as 



between the ring and the reservoir a can in general be 
expressed as 

4 C = G v dc v + G a dc n + G^nl + GT x vn , (9) 



which satisfies J'jf 



- J^ c , as is expected from the con- 
servation of the charge. In the above equation we define 
four different transport coefficients. The coefficient G\ c 
is the usual linear dc conductance, G a dc is the coefficient 
relating the dc current with the pumping frequency in 
the adiabatic pumping regime, G™ is the non-adiabatic 
transport coefficient and G™ 1 is a coefficient that quan- 
tifies the effect of mixing between the dc bias and the 
pumping to generate the dc current. The corresponding 
expressions are 



E En^(M)r a W|5j [ tl >,m)I 2 , 

- £ Znr fi (ji)T a (ji) d u \g?,(n,uj)\ a 

P=L,R n 

G dc = E Er^)r„M w)| 2 



Ga 
dc 



Gmix 
dc 



LO—fJ, 

(10) 



n2Re 



where the matrix elements Q\ l (n, uj) with k = 0, 1 cor- 
respond to a low- frequency expansion of the Green's func- 
tion of the form 



g(n,u) ~g°(n,w) + ^ g 1 (n,w). 



(11) 



In the next sections, we will analyze the dependence of 
the non-vanishing pumping coefficients on the parame- 
ters and symmetries of the setup. On the other hand, we 
will also discuss the symmetry properties of the current 
as a function of $^ c , j = 1,2. 



H r (t)=H +V(t), 



(6) 



where Ho and V(t) = En^o e~ ma °*V„, are matrices with 
elements defined by the spacial coordinates of the rings. 
We then formulate the Dyson equation as follows 

G R {t,uj) =<5°(w) + 

J2e- innot G R (t,u} + nn )V n G°(ij), (7) 

n=£0 



C. Evaluating the retarded Green's function 

In this subsection we present the different strategies 
that we will follow to evaluate the retarded Green's func- 
tion entering the expression for the current in different 
relevant limits. 



Small ac amplitudes 



where G R (t,u) denotes the matrix with elements 
G R ll {t,ui), being 1,1' spacial coordinates of the rings, 
while the stationary retarded Green function is also ex- 
pressed as a matrix 



G°(w) = [ui-ffo-EH]" 1 , 



(8) 



with Z ltl ,(u})=J2 a 5 lt t,5 lt i ct J{dLj')/(2w)r a (cj')/(Lj-Lj' + 
t0+). 

For a very small voltage difference V and low driving 
frequency Oq, the dc current flowing through the contact 



The solution of the Dyson equation leads to the 
exact Green's function. For weak amplitudes 3>f c and 
$2 C of the ac component of the magnetic flux and arbi- 
trary frequency J7o, it is possible to solve that equation 
perturbatively. It is convenient to start expanding the 
Hamiltonian H r (t) in powers of <&1 C , $% c . By keeping 
terms up to the first order in these parameters, we get 

2 N 

H r (t) ~H +J2J2 [v^dHjCl+ij + H.c] , (12) 
j=x i=i 
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with 



and 



with 



V {j \t) = i Wj e-^T(j>f cos(fi * + Sj), (13) 



with <fi dc ' ac = $" c ' ac /L, j = 1,2. The corresponding 
perturbative solution of Eq. ([7]) reads 



2 iV 



-id> da f 



^dc,ac 



(14) 



G R (t,iu) ~G°(w) + e ~ mn0t 

ri=-l,n/0 

xG°{u + nn )V n G (uj), (15) 
where V™ = V">W + V n '( 2 ), with matrix elements 



6i> 



i',i+i e 



J = 1,2, 
(16) 

being 1,1' sites of the ring j = 1,2, while G°(w) is given by 
Eq. ([5]). This procedure can be sistematicaly extended 
to consider higher order solutions in $° c , $2 C . 

2. Single ring weakly coupled to reservoirs 

For the case of a single ring, a possible route to calcu- 
late the retarded Green's functions, alternative to solving 
([7]) consists in starting from the limit where the ring is 
completely uncoupled from the reservoirs. The corre- 
sponding Hamiltonian can be recasted as 



k Ck, 



(17) 



e k {4> dc ) = -2wcos(k + (j) dc ), 
v k {4> dc ) = 2wsm(k + <p dc ). 



(21) 



The Green's function including the coupling to the leads 
is the solution of the following Dyson's equation 



G(m,uj) — g(m,uj) + ^£(m — n,u + nQ, ) 

n 

xS(u + nflo)g(n, w), 



(22) 



where the matrix g(n,u>) has matrix elements giji(n,uj). 
This equation can be exactly solved or can be used to ob- 
tain perturbative solutions in the coupling to the reser- 
voirs and the strength of the potential profile. 



D. Low- frequency expansion 

For low frequencies a solution exact up to O(£lo) can 
be obtained by expanding Eq. ([7]) as follows: 



G R {t,u) ~ G (0) H + G R {t,Lo)V(t)G (0 \co) + 
z^G*(^)^G(°>H. 



We define the frozen Green's function 



G f (t,iu)= G^iu})- 1 - V(t) 



(23) 



(24) 



in terms of which the exact solution of the Dyson equa- 
tion at 0(f2o) reads 



where c k = l/v^/V^i^i e lkl °U with k — Ivn-njN, with 
-N/2 < m < N/2 and e k (0(f)) = -2w cos (fc + </>(*)). 
The exact retarded Green function for this problem is 



g R (t, f') = -iQ(t - t') cxp{-i / dhsk (0(ii))}.(l 



In the limit of small $ ac , keeping terms up to the first 
order in this parameter, we can express 

g R (t,t') = £ e-^t f ^ e -*-(t-O flfc(n(a>)) (19) 



n=-l 



2vr 



being 



g k {0,uj) = g° k (uj) 



1 



g k {±l,u) = ± 



uj - e k (<p^) + ir\ 

r C V k (t dC ) r o 



2ft 



[g k ( U )-g k (u>±n )] ,(20) 



G (1) (f , u) = & (t, u) + id u G f (f , tv)^jp-G f (f , u). (25) 

As we will discuss in the next section, this approach is 
useful to evaluate the transport coefficients in terms of 
the frozen Green's function within the adiabatic regime. 
This will allow us to analyze the symmetry properties of 
the current as a function of $^ c , and for different config- 
urations of the attached reservoirs. 



III. CONDITIONS FOR ADIABATIC PUMPING 

Our first step is to analyze the conditions to get adia- 
batic pumping, which implies a non- vanishing coefficient 
G% c . Using the Floquet- Fourier representation of Eq.© 
in Eq. ([25]) . we can identify the first term of the expan- 
sion (Ull) 



Jo T 



(26) 
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where r = 27r/£7o is the oscillation period of the magnetic 
flux $° c (t). Replacing this expression in the adiabatic 
coefficient of Eq . (flOf . we obtain 



G a dc = E r^( M )r Q ( M )x 

0=L,R 



(27) 



where from Eq. (|24|) we can calculate the derivative of 
the frozen Green function 



d t &{t,w) = G f {t,Lu)d t V{t)G f {t,Lu). 



(28) 



It was shown in Ref. |22j that at least two parameters 
are necessary to have adiabatic pumping in ac driven 
systems. In what follows we show that, in the present 
context, this implies at least two different rings driven 
by two different magnetic fluxes. 



A. Single ring 

The first question that arrises is about the possibil- 
ity of implementing a driving with a magnetic flux char- 
acterized by two or more parameters in a single ring. 
Such a possibility would correspond to a single ring 
threaded by a flux containing several harmonics, of the 
form $(t) = $ dc +£m=i $ ac >( m ) cos{mn t+S m ). In such 
a case, we can express 



H r (t) = e m w + e~ i<Kt W, 



(29) 



with (f>{t) = ${t)/L and W^ v = -wSfj+i, with / = 1, N 
and N + 1 = 1. 

In the appendix A we show that the frozen Green's 
function in this case has the following structure 

G f (t,u)= e mm G {n) {uj). (30) 



with 
A (1) 

A (2) 

with 



^ I>(/i)r a (/i)Im A 1 (n)i a , h A 2 (fi) 



P=L,R 



£ Wr^Im [AsO^^A^ J(,33) 

P=L,R 



(34) 



where Jj is the current operator along the circumfer- 
ence of the ring j, defined in (|B2p . These equations 
allow us to identify two conditions for a non-vanishing 
adiabatic current. In fact, after performing the calcu- 
lations of Eq. (1321) explicitly, it is easy to verify that 
G% c oc 4>1 c (t>2 C sin(<$i — <5 2 ). Thus, the first condition is 
a phase difference for the ac fluxes, Si — 62 =/= nir. The 
other condition is determined by requesting non vanish- 
ing matrix elements Aj(/x); a) ;_, which implies dc mag- 
netic fluxes threading both rings satisfying 7^ nir. 
These two conditions are in agreement with the results 
of Ref. 3I- In Fig. [2] we show that these conditions also 
hold for arbitrary (not necessarily small) amplitudes of 
the ac components. We observe that as the amplitude of 
the ac components increase, G a dc as a function of S± — 82 
shows a much more complex structure than the simple 
sinusoidal law predicted by the small amplitude result of 
Eq. (f3"2"j). However, the condition of a vanishing value of 
this coefficient when the phase difference coincides with 
an integer value of tt is verified in all the cases. As a 
function of the magnetic flux, it is also observed that for 
increasing ac amplitudes, there are sign changes in the 
behavior of the pumped current, but this current van- 



ishes for $ 



dc 



nir. Further details of the behavior of 



Gj c as a function of will be analyzed in the next 
section. 



The integral in time entering (1271) is proportional to 



d$(t) 
' — — 1 
dt 



i{n - m)Ht) G^iM\ G t\S^\\ (31) 



which vanishes for any periodic $(t), implying a van- 
ishing adiabatic coefficient G^ c . We, therefore, conclude 
that it is not possible to generate an adiabatic pumped 
current by applying pure harmonic magnetic fluxes in a 
single ring. 



B. Two rings 

Using the frozen Green's function for small ac ampli- 
tudes of the appendix B in Eq. (JBTJ) leads to the following 
adiabatic coefficient 



2vr Jo dt 2tt J 



wry 

dt ' 

(32) 



IV. DEPENDENCE OF THE PUMPED 
CURRENT ON THE STATIC MAGNETIC FLUX 

A. Single ring 

As discussed in the previous section, it is not possible 
to have in this case a pumped current within the adi- 
abatic regime. The non-adiabatic transport coefficient 
G™ is, however, non- vanishing and we now turn to study 
its behavior as a function of the dc magnetic flux. In par- 
ticular, we are interested in analyzing if the non-adiabatic 
pumped current has a defined parity as a function of the 
dc magnetic flux, as in the case where Onsager-Casimir 
relations are valid. 

In the limit of a weak coupling between the ring and 
the reservoirs and for small amplitudes of the ac fluxes, 
it is possible to find an analytical expression for G™- 
Evaluating the dc current J^ c at the lowest order in the 
couplings |w Q | 2 , corresponds to considering S(w) — > in 



6 



cient (flOl) results 
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FIG. 2. Top panel: G^ c as a function of the dc component of 
the magnetic flux $ dc = $f c = $a c , for two rings with TV = 10 
sites connected at h = 6 and h = 1 with w c — 1. The leads 
are symmetrically coupled to the ring with w\ — 0.5 and 



0.2. The chemical potentials are the same for the two 
reservoirs, fi = 0.18. Energies are expressed in units of the 
hopping matrix element wi = W2 = w. The ac flux of the left 
ring is fixed at the value 0J C = 0.2, while 4>2 C is varied and the 
phases of the ac fluxes are <5i = and 82 = vr/2. Solid line, 
triangles, dot-dashed line, dashed line and circles correspond, 
respectively, to </>g c = 0.45, 0.4, 0.35, 0.3, 0.2. Bottom panel: 
G a da as a function of the phase difference of the ac fluxes 82 — 
Si = 8. The ac components of the fluxes are <j)i c = 0.2 and 
4>2 C = 0.5, and the dc fluxes of the rings are equal. Dashed 
and solid lines correspond to <J> dc = 1.4, 3.8. 



I, which implies Q{rn,uS) ~ g(m,uS). For small fioi 
performing the expansion of Eq. in these Green's 

functions cast 



e»^(±i, w ) = -^, t (^)fe-^-w 



duj 

2„0 



^( ±i . w )==fSe^)^'-- , ' , < 35 ) 



Replacing these expressions in the non-adiabatic coeffi- 



2Re 



duj duj 2 



where II and Ir are the sites of the ring where the left 
and right reservoirs are connected. The dependence of 
this coefficient on the dc magnetic flux is through the 
energies Ek(4> dc ) and the currents Vk(4> dc ) given in Eq. 
(|21jl. In the case of reservoirs symmetrically coupled to 
the ring, we have Ir — Il = N/2, which corresponds to 



i(k-k')(l L ~l R ) 



The sums in (1361) are, thus in 



variant under the change k 

e- k {-^p dc ) = e k ((j) dc ) and V- 



-k and k! —> —k'. Since 

(j) dc ) = — Vk ^(j) dc ^ we can 



conclude that the non-adiabatic coefficient is an even 
function of the dc magnetic flux § dc when the reservoirs 
are symmetrically connected. However, when the reser- 
voirs are connected at arbitrary positions the phase fac- 
tor e l ( k ~ k )('^ _/ «) is no longer invariant under inversions 
of k. Then, the coefficient does not have a defined 
symmetry as a function of $ dc . 

In Fig. [3] we show the behavior of the non-adiabatic 
pumped current as a function of $ dc for a ring threaded 
by a flux with also an ac component < I> ac and reservoirs at- 
tached at symmetric and asymmetric positions obtained 
by numerically solving the Dyson equation ([7]) in the 
limit of a small <f> ac . In order to have a non- vanishing 
pumped current, it is necessary to break the spacial in- 
version symmetry^ which in our case is accomplished 
by connecting the reservoirs with different tunneling am- 
plitudes wl 7^ wr. These results show that the above 
conclusion on the behavior of the pumped current as a 
function of <5> rfc is also valid for arbitrary couplings be- 
tween the ring and the reservoirs. In fact, the pumped 
current is an even function of $dc when the reservoirs are 
coupled at symmetrical positions along the ring, and it 
has no particular symmetry when the coupling is asym- 
metrical. 

Another remarkable feature that is observed in some 
cases with asymmetric coupling to the reservoirs (see for 
instance the plot in circles of Fig. [3]) is the fact that, 
as a function of the driving frequency Oo; the dc cur- 
rent changes from a paramagnetic-like behavior, charac- 
terized by a vanishing magnitude at <& dc = at low fio 
to a diamagnetic-like behavior, characterized by a siz- 
able amplitude at <f> dc = 0. We have verified that such a 
change in the behavior takes place for a driving frequency 
CIq ~ 0.2, which corresponds to resonance with the mean 
level spacing of the ring and it is then likely to be caused 
by interference effects introduced by the ac driving. 



B. Two rings 

We showed in the previous section that, in the case 
of two rings it is possible to generate a finite pumped 
current within the adiabatic regime. The corresponding 
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FIG. 3. The pumped current as a function of the dc com- 
ponent of the flux $ dc for a single ring with N = 20 sites, 
when flo = 0.2. The leads are symmetrically coupled to the 
ring with w\, = 0.5 and w R = 0.2. The chemical potentials 
are the same for the two reservoirs, = Hr = 0.18. En- 
ergies are expressed in units of the hopping matrix element 
w. The ac component of the magnetic flux is <f) ac — 0.01 and 
5 = 0. Solid line corresponds to wires coupled at symmetrical 
positions of the ring, while dashed line and circles correspond 
to asymmetric coupling with Ir — II = 4 and Ir — Il = 7, 
respectively. 



transport coefficient is given in Eq. Q32p and we now 
focus on its symmetry properties as a function of the 
dc components of the fluxes <£>^ c . The dependence of 
G a dc on Qj c is enclosed in the coefficients through 
the matrix elements Aj(/j,)i a ^ p . This matrix is defined 
in Eq. (|34[) . The kinetic energy operator of each ring, 

Wje 1 ^" + Wje~ l ^ c entering the Green's function G°(w) 
is an even function under the transformation Sj : </>j c — > 
—(j)j C , ( — h j)i corresponding to a simultaneous 

inversion of the flux (f>j C and a spacial inversion along the 

circumference of the ring, while the current operator Jj 
is odd under such transformation. Notice, however, that 
the contacts between the rings as well as the contacts 
between the rings and the reservoirs also enter G a (ui). 
As a consequence, Sj are symmetries of the full setup 
provided that the couplings between the rings and/or 
between rings and reservoirs do not break them. 

Interestingly, for symmetrically connected rings and 
reservoirs, the transport coefficient G a dc and the corre- 
sponding adiabatic current are odd functions of each of 
the dc fluxes, considered independently one another, as 
illustrated in the top panel of Fig. @J while they are even 
functions of the dc magnetic flux when it is simultane- 
ously varied in the two rings, as shown in the bottom 
panel of Fig. 2J For arbitrary couplings between the 
rings and reservoirs, the pumped current does not have 
any particular symmetry as a function of the dc fluxes, 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 4. The adiabatic coefficient G^ c for two rings driven 
by harmonic fluxes and symmetrically connected to the reser- 
voirs, as functions of the dc components of the fluxes. In 
both panels, the phases of the ac fluxes are Si — and 
5-2 = 7r/2. Top panel: the dc flux of the right ring is kept 
fixed and the dc flux of the left ring $f c is varied. The ac 
fluxes are <f>i c = (j>2 C = 0.2. Dashed line, circles, solid line 
and triangles correspond, respectively, to $2 C = 0, 0.6, 3, 5.3. 
Bottom panel: the two dc fluxes are simultaneously changed 
$f= = <I>^ C = $ dc . Solid line, circles, and dashed line corre- 
spond to 4>T — 4>2 C = 0.4, 0.3, 0.2. Other parameters are the 
same as in Fig[2] 



meaning that Onsager symmetry is not expected to be 
observed in this general case. This is illustrated in Fig. [S] 
where the adiabatic transport coefficient is shown for the 
case of two rings asymmetrically coupled to reservoirs. 



V. DC CURRENT WITH AC FLUXES AND 
BIAS VOLTAGE 

We complete the analysis of the symmetry properties 
of the dc current as a function of the dc magnetic flux by 
adding the effect of a bias voltage applied as a chemical 
potential difference at the two reservoirs. 



8 




l2 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



<1> 



FIG. 5. The adiabatic coefficient G% c for two rings driven by 
harmonic fluxes asymmetrically connected to the reservoirs, 
as functions of the dc components of the fluxes. The two dc 
fluxes are simultaneously changed <J>f c = <E>2 C = $ dc . The ac 
components of the fluxes are 4>t c = 0.2 and ffif = 0.4, and 
the phases are 5i = and 82 = n/2. Circles corresponds to 
wires coupled at II = 10 and Ir — 2, dashed line corresponds 
to II — 5 and Ir = 9, while solid line corresponds to II = 1 
and Ir — 2. Other parameters parameters are the same as in 
FigH 



A. Single ring 

In this case, we have shown in Section III that the 
adiabatic coefficient is G a dc — 0. A similar analysis cast 
G™^ = 0, which means that for small V and f^o, the 
voltage and the pumping contribute independently to the 
dc current. It is well known that the linear stationary 
conductance obeys Onsager-Casimir relations, irrespec- 
tively the details of the contacts to the reservoir. Thus, 
from the analysis of the previous section we conclude that 
the full dc current of the driven ring presents Onsager- 
Casimir symmetry only for the case of reservoirs that are 
symmetrically connected. 

In Fig. [6l we show the behavior of the total dc current 
Jdc = G dc v + J pump , as a function of the dc magnetic 
flux, for a ring with symmetrically connected reservoirs 
under the combined effect of a small bias voltage and a 
magnetic flux oscillating with a small amplitude $ ac and 
several frequencies Qq beyond the adiabatic regime. For 
the lowest frequencies jv um v ^ G™&o, a s discussed in 
section III. In the upper panel we show the linear dc con- 
ductance G^ c , which is an even function of $ dc , in the 
middle panel the total current J dc /V. In the latter case 
we divide the current by the voltage in order to show a 
quantity having the same units as the usual linear con- 
ductance G dc shown in the top panel. For the symmet- 
ric connection, the pumped current JP um P for vanishing 
bias voltage V = is an even function of for all the 
frequencies considered (see the botton panel of Fig. [B]) . 



Therefore, the total current J is also an even function 
of $ dc . The corresponding behavior for asymmetric con- 
nections of the reservoirs is shown in Fig. [7J In this case, 
although the linear conductance G dc shown in the top 
panel is an even function of $ , the pumped current, 
shown in the bottom panel, does not have a well defined 
symmetry. Thus, the total dc current shown in the mid- 
dle panel of Fig. [TJdoes not have a well defined symmetry 
as a function of $ dc either. 

Notice that the quantity J dc /V shown in the middle 
panels of Figs. [5] and [7] can be significantly larger than 
the conductance quantum expected for a single channel 
system, like the ones shown in the top panels. This be- 
havior is precisely a consequence of the fact that the dc 
current is not only induced by the bias voltage but also 
contains a non-vanishing component due to the ac driv- 
ing. 
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FIG. 6. Upper panel G dc as a function of the dc component of 
the magnetic flux $ tic , for a single ring with N = 20 sites. The 
chemical potentials are /j,r = fj,+V and /j,r = \i, with /i = 0.18 
and V = 0.001. The reservoirs are symmetrically attached to 
the ring, and coupling parameters are the same as in Fig. [3] 
In the middle and lower panels dashed line, circles, solid line 
correspond, respectively, to Qo = 0.08, 0.14, 0.2. Lower panel: 
pumped dc current for = Hr = A 4 - The dc current resulting 
from the combined effect of ac driving and dc voltage bias is 
shown in the middle panel divided by the voltage bias. Other 
parameters are the same as in Fig. [3] 



9 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 7. The same as Fig. [6] for reservoirs asymmetrically 
attached at II = 1 and Ir = 10. 

B. Two rings 

In the case of two rings under the combined effect of ac 
driving and dc bias, the full dc current for small V and Qq 
can be expressed as J dc = G% C V + Q, Q G a dc + Q. Q VG 7 £ x . In 
terms of the frozen Green's function the mixed transport 
coefficient can be writen as 

GT X = - £ Wr a ( M ) 

P=L,R 

^ a (37) 
Jo £ ^Im (G{ ath (t, u)d t G{l h (t, u)) . 

In the limit of small amplitudes a similar treatment 
and analysis to the one performed in section III leads 
us to the conclusion that this transport coefficient has 
the same symmetry properties as a function of ^ c as 
G a dc . Namely, for a symmetrical configuration of rings 
and reservoirs, this coefficient is odd under inversion of 
one of the fluxes and even under the simultaneous inver- 
sion of the two fluxes. On the other hand, G^ c is in this 
case an even function under the inversion of any of the 
fluxes &j C , and also under the simultaneous inversion of 
the two fluxes. Therefore, in the presence of a bias volt- 
age and in a setup symmetrically connected, we expect a 
total dc current with no particular symmetry when a sin- 
gle flux is inverted, while we expect a dc current which is 



even under the simultaneous inversion of the two fluxes. 
This is illustrated in Fig. [5] where we consider ac fluxes 
with the same amplitudes driving both rings. We fix the 
dc component of the flux threading the right ring and 
analyze the dc current as the flux threading the left ring 
changes. We show in the top panel of the figure the dc 
current resulting from the bias voltage divided by the 
voltage V, G v dc + G^ X VL . The middle panel shows the 
mixed coefficient G™ x and in the bottom panel the total 
current divided by the voltage J dc /V. The correspond- 
ing adiabatic coefficient has been shown in Fig. 2) In 
agreement with the analysis for small ac amplitudes, no 
particular symmetry of the total current as a function of 
the flux <&f c is observed, and this behavior is not consis- 
tent with Onsager-Casimir relations. In Fig. [H] we show 
results for the same setup considered in Fig. |5]but we now 
change the two fluxes simultaneously, <&f c = $2 = $ dc . 
We see that in this case, all the transport coefficients are 
even functions of 4> dc (recall that G a dc is shown in Fig. @|. 
The behavior of the total current shown in the bottom 
panel of the Fig. [3] is in the present case consistent with 
Onsager-Casimir relations. 

In the case of two rings asymmetrically connected, G dc 
is an even function of $^ c , j — 1,2, but the adiabatic 
component does not have any defined parity (see Fig. [5]) 
and we have verified that this is also the case of the mixed 
component. Thus, the full dc current does not have in 
this case any particular symmetry as a function of dc 
component of the magnetic flux. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the transport properties of one and 
two rings threaded by magnetic fluxes with dc and ac 
components with and without dc bias voltage applied at 
the terminals. We have developed different theoretical 
strategies to solve the problem in different limits and we 
have defined the relevant coefficients to characterize the 
transport in the limit of small bias voltages and low driv- 
ing frequencies. 

We have shown that it is not possible to generate adi- 
abatic pumping in a single ring by driving with pure ac 
magnetic fluxes, even for magnetic fluxes containing sev- 
eral harmonics oscillating with phase-lags. This is in con- 
trast to the behavior of a single ring threaded by a mag- 
netic flux that changes linearly in time. In that case, 
the induced constant electric field is proportional to the 
driving frequency and the induced dc current is, thus, 
proportional to the frequency even when it corresponds 
to a single-parameter driving.* 9 In the case of two rings 
we have shown that adiabatic pumping is possible pro- 
vided that the two magnetic fluxed have ac components 
oscillating with a phase lag and also have a finite dc com- 
ponents, which is in agreement with previous results^ 

Finally, we have analyzed the behavior of the dc cur- 
rents in setups with one and two rings as the dc magnetic 
flux is varied with and without applied dc voltage. For 
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FIG. 8. Upper panel: The dc current resulting from the dc 
bias voltage divided by the voltage V, Gj c + G™ x Q,q, for 
two rings symmetrically connected. The dc component of 
the magnetic flux of the left ring $f c is varied and the one 
threading the second ring is kept fixed $f = 0.6. Middle 
panel: The transport coefficient G mrx . Lower panel: the total 
dc current resulting from the ac driving combined with the 
effect of the dc bias voltage divided by the voltage V, when 
fio = 0.02. The dc bias voltage is V = 0.001, and the ac 
components of the fluxes are (j>i c = 4>2 C = 0.2. The phases of 
the ac fluxes are Si = and 82 ~ n/2. Other parameters are 
the same as in Fig. [21 



the case of a single ring we found that the non-adiabatic 
pumped current does not have in general any particular 
symmetry as a function of the dc flux, indicating that 
the Onsager-Casimir relations are not valid in this sys- 
tem. An arbitrary small coupling to the reservoirs in 
the driven ring is enough to break Onsager symmetry, 
except in the case where the connection is at perfectly 
symmetric positions under spacial inversion symmetry. 
We found a similar behavior in the case of two rings and 
pumping within the adiabatic regime. For the particu- 
lar case of perfectly symmetric coupling to the reservoirs, 
the pumped current in this setup is even as a function of 
the dc magnetic flux, when the same dc flux threads the 
two rings, while it is an antisymmetric function under 
changes of the dc flux of only one ring while keeping the 
constant the dc flux of the other one. The fact that in 
the non-adiabatic as well as in the adiabatic regimes the 



FIG. 9. The same as Fig. [8] for simultaneous variations of the 
dc magnetic fluxes of the two rings $f c = <E>2 C = $dc- 



Onsager-Casimir relations are not expected to be valid 
in general can be related to the fact that pumping is at 
least a second order process in the driving amplitudes as 
explicitly shown by Eq. (|36p for a single ring and Eq. 
([32]) for two rings. Although in the adiabatic regime the 
dc current is linear in the pumping frequency it is non- 
linear in the driving field. In this sense, the situation 
resembles the case of the dc driving where Onsager sym- 
metry is broken beyond linear response in the applied 
voltage. However in that case, the explanation of this 
breakdown resort to the effect of the interactions^ while 
in the present case, it is a purely dynamical effect, which 
takes place in a non-interacting system. When a dc bias 
voltage is, in addition applied, the dc current resulting 
as a combination of driving with the ac flux and with 
the dc voltage is an even function of the dc magnetic 
flux only in the case where the setup is symmetrically 
connected, not showing particular symmetries in other 
configurations, in agreement with the breakdown of the 
Onsager-Casimir relations in the pumped component of 
the current. 

To finalize, we would like to mention other inter- 
esting features of the transport properties of these se- 
tups. In particular, in the case of pumping in a single 
driven ring, we would like to stress the change from the 
"paramagnetic" -like behavior, characterized by a vanish- 
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ing dc current for vanishing dc flux, to the " diamagnetic" - 
like behavior, characterized by a finite current for zero dc 
flux , as the driving frequency increases. This feature is 
akin to what has been experimentally observed in the 
behavior of the persistent currents in rings threaded by 
fluxes with ac components.— 
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Appendix A: Frozen Green's function for a single 
ring 

We consider a single ring threaded by a magnetic flux 
with an arbitrary number of harmonic components de- 
scribed by the Hamiltonian (|29|) and connected to reser- 
voirs. The frozen Green's function can be expressed in 
terms of the T-matrix as 

Gf(t,cj) = gf(u,)+gf(uj)ff(t,u,)gf(u,), (Al) 



being gf(u) 



A - E(oj) 



, and 



f f (t,Lu) = H r {t)+H r (t)g f (uj)H r (t) + 



Appendix B: Frozen Green's function for two rings. 
Fluxes with small ac amplitudes 

For small amplitudes of the ac components of the 
fluxes, we can consider the Hamiltonian (|12|) to find that 
the frozen Green's function at the first order in <j>j (t) 
reads 

2 

d f (t,u) =G (u) + Y,h( t )G (u)jjG°{cj), (Bl) 
i=i 

being GVr 1 = - H - £, with H defined in (|T3|) 
while 



(B2) 



(A2) 



defines a matrix associated to the persistent current op- 
erator of the ring j. 
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